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Abstract 

Quantum electrodynamics has been the first theory to emerge from the ideas of regulariza- 
tion and renormalization, and the coupling of the fermions to the virtual excitations of the 
electromagnetic field. Today, bound-state quantum electrodynamics provides us with ac- 
curate theoretical predictions for the transition energies relevant to simple atomic systems, 
and steady theoretical progress relies on advances in calculational techniques, as well as 
numerical algorithms. In this brief review, we discuss one particular aspect connected with 
the recent progress: the evaluation of relativistic corrections to the one-loop bound-state 
self-energy in a hydrogenlike ion of low nuclear charge number, for excited non-S states, 
up to the order of a (Za) e in units of the electron mass. A few details of calculations for- 
merly reported in the literature are discussed, and results for 6F, 7F, 6G and 7G states are 
given. 

PACS Nos.: 31.30.Jv, 12.20.Ds, ll.10.St 
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1 Introduction 

The bound-state self-energy of an electron in a hydrogenlike atom has been one of the key theo- 
retical and experimental cornerstones in the development of modern field theory. During the last 
decade, the techniques available for calculations have dramatically advanced, and the accuracy 
of today's predictions is several orders of magnitude better than in the early days of quantum 
electrodynamics (QED). The advances would be impossible without the combined theoretical 
efforts toward the description of the hydrogen spectrum, which have involved generations of 
physicists. The aim of the current brief review is threefold: first, to give an overview of some 
recent advances in Lamb shift calculations and self-energy calculations in particular, second, 
to describe a few details of recent calculations, for higher excited non-S states, which may be 
useful in an independent recalculation of the higher-order binding corrections, and third, to sup- 
plement previously available data for higher-order corrections by results for 6F, 7F, 6G and 7G 
states. 

In "usual" calculations of virtual loops, we are used to the association: "the number of loops 
equals the power of the coupling constant." For hydrogenic bound states, the situation is dif- 
ferent, and the expansion of a bound-state energy, including QED effects, is actually a double 
expansion in terms of the QED a fa 1/137.036 and the electron-nucleus coupling parameter 
Za, where Z is the nuclear charge number. 
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Specifically, the leading one-loop energy shifts (due to self-energy and vacuum polarization) 
in hydrogenlike systems are of order a (Z a) 4 mc 2 , where m is the electron mass and c is the 
speed of light. The complete correction in this order was obtained first in the years 1948 and 
1949 yj,l2l|3l|4j]. The leading two-loop energy shifts are of order a 2 (Z a) A mc 2 \ the paradigm is 
that a counts the number of loops, whereas Za is a measure of the relativistic corrections that 
enter the loop corrections and are typical of the bound-state problem. 

When scanning the literature, one should be aware that sometimes, the QED energy shifts are 
expressed in terms of atomic units, for which the fundamental energy scale is that of the Hartree, 
which equals a 2 mc 2 , and therefore the leading QED shift is of order Z 4 a 3 in atomic units. 

The current paper will be concerned mainly with the calculation of the correction of order 
a (Z a) 6 , which is effectively a relativistic correction to the one-loop correction (sic!) due to 
the binding of the electron to the nucleus (hence the name "binding corrections"). We thus have 
to consider the relativistic atomic physics aspect of the problem in addition to the usual loop 
calculation. 

The separation of the problem into high- and low-energy virtualphoton contributions has been 
discussed (e.g.) in §123 of Ref. [5]. In the context of NRQED [atZllElSl], the general paradigm is 
to map the high-energy effects onto effective operators. By contrast, the low-energy effects are 
mapped onto the transverse degrees of freedom of the electromagnetic field and are integrated 
up to some cutoff scale. The systematic expansion of the Wilson coefficients multiplying the 
effective operators is rather nontrivial. In the current context, we start from the fully relativistic 
expression of the self-energy, and do appropriate expansions for both domains. For the problem 
at hand, this conceptually simpler approach has certain pragmatic advantages. 

One might well ask why we would need information on the Lamb shift of higher excited states in 
the first place, especially because the most accurately measured hydrogenic transition involves 
the IS and the metastable 2S level [10] . The answer is that one transition is not enough to 
determine fundamental constants like the Rydberg to sufficient accuracy: one needs at least two 
transitions. We would like to cite from Ref. Ill ill , a system of equations that illustrates this fact 
in a particularly clear way. One combines two of the most accurately measured frequencies 
(/is-2S ifioll and a /2S-8D 1121 transition), where we can either use 8D 3 / 2 and 8D 5 / 2 for the 



8D level, for a determination of the Rydberg. Defining eo as a dimensionless relativistic Dirac 
energy, we can schematically establish the following equations, 

/is-as =#00 c {e D (2S) - e D (lS)} + C 2S - £ 1S , (la) 

/2S-8D =Roo c {e D (8D) - e D (2S)} + £ 8D - C 2S • (lb) 

By L, we refer to the Lamb shift [see Eq. (67) of Ref. iflill for an often-used definition of this 
effect]. We can eliminate L 2S using the combination 

Lis + W 2l 

^2S = g lA) 

where 

W 21 = 8L 2S - Lis (3) 

is independent of the nuclear-size correction (proton radius) and can therefore be determined 
theoretically to high accuracy. We can now solve the system of equations <[TJ for two unknowns: 
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Lis an d -Rqo, provided we have theoretical values for the weighted combination W 2 \ and L 8D . 
The evaluation of an important contribution to the latter will be one of the issues discussed in 
the current brief review. 



A generalization of this rather elementary exercise involves a least-squares analysis [14]. 
Roughly speaking, one tries to find the "best fit," in the sense of least squares, to a set of 
experimental results representing the most accurately measured transitions in hydrogen and 
deuterium. For the 2002 adjustment [15], a total of 23 transitions in hydrogen and deuterium 
have been used as input data. The least-squares analysis includes covariances among the the- 
oretical uncertainties. One thus obtains theoretical values, including standard uncertainties, for 
the proton radius and the Rydberg, and one can make optimal predictions for other transition 
frequencies not used in the adjustment of the constants. Some of the predicted frequencies are 



currently more accurate [14] than the predictions for the anomalous magnetic moment of the 
electron. 

There is another aspect which is relevant to calculations for higher excited states: the extrapo- 
lation of the data to higher quantum numbers, for given manifolds of states. For example, we 
can extrapolate to higher principal quantum numbers n for given angular momenta (orbital and 
spin). In an apparently not widely known paper (Ref. lfl^1 \ the one-loop Bethe logarithm has 
been evaluated, in terms of a (slowly convergent) integral representation, in the limit of n — > oo. 
From this point onwards, the extrapolation to higher principal quantum numbers of the one-loop 
nonrelativistic Bethe logarithm, based on a few known low-n results, becomes an interpolation. 
However, this calculation uses a few particular properties of the nonrelativistic Bethe logarithm 
[order a(Za) 4 ], which are not applicable to the much more involved correction of order a (Za) 6 
(relativistic Bethe logarithm). Nevertheless, it has recently been observed that a large number 
of relativistic effects and loop corrections can actually be expressed as an inverse power series 



in the principal quantum number [ 17, 18]. This is not irrelevant because it gives confidence that 
such fits, applied to other corrections like the relativistic Bethe logarithm provide us with an 
accurate representation of the corrections for large n, where the complex structure of the wave 
function leads to a prohibitively large number of terms in intermediate steps of the calculation, 
which makes an evaluation impossible even if today's computer algebra systems f 1911 are used. 

This brief review follows two investigations in which previous treatments of the one- 

loop problem were extended to higher quantum numbers, and some general expressions rel- 
evant to selected manifolds of bound states were presented. Here, we supplement the results 
obtained previously by explicit numerical calculations for 6F, 7F, 6G and 7G states. Further- 
more, in Sec. |2j we discuss some general ideas relevant to the calculation, before dwelling on 
the low-energy part in Sec. |3] This part, mediated by ultrasoft virtual photons, is one of the 
key elements in the treatment of higher excited non-S states. Throughout the paper, we restrict 
the discussion to states with nonvanishing orbital angular momentum. Numerical results are 
presented in Sec.HJ and conclusions are drawn in Sec.|5] 



2 General Ideas 

We first take the opportunity to refer the reader to the two most comprehensive treatments of 
the one-loop self-energy problem for a bound state, which have appeared in the literature up to 
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now and which we are aware of. These are Refs. [21] and 12211 . Here, we briefly mention a few 
particular aspects of interest which may be useful in finding a general access to the calculations, 
also with regard to a possible independent verification. We fix the units so that h = c = e = 1, 
which implies that e 2 = Area. The electron charge is e, not — e. 

In Feynman gauge, the one-loop bound- state self-energy reads 



E 



i e 



d 4 k 



1 



r 



7/i 



ip j — (ij) \5m\ ip^ 



(4) 



where D fll/ (k) is the photon propagator, and Cp is the Feynman contour. Lorentz covariance is 
explicitly broken by the Coulomb potential V = —Za/r. The quantity 5m is the mass counter 
term, and m is the electron mass. 

We divide the problem of calculation into two contributions, which correspond to two different 
photon energy regions, and two different parts of the photon energy integration contour Cp. 
The photon energy separation parameter e takes the role of an infrared cutoff for the high- 
energy part, and it acts as an ultraviolet cutoff for the low-energy part. The low-energy part may 
be expressed as 

E L = aJ2(Za) n f n f \ 9 j . (5) 

n>4 ^ 



(Za) 2 m 



and the high-energy part E H reads 



E H = aJ2(ZaTg n (£) . 



(6) 



n>4 



The key is to first expand both contributions in Za for small Za, then in e, for small e. By per- 
forming the expansions in that sequence, we automatically assume that (Za) 2 m <C e. Counter- 
intuitively, this is equivalent to performing an expansion for large e, after scaling the powers 
of Za out of the calculation, as will be demonstrated below in Sec. 13.31 It is also important 
to realize that the actual numerical value of e is arbitrary. All that matters is the expansion in 
systematic expansion in e, including the logarithmic terms. In Refs. [21] and H22I1 . it has been 
stressed that different gauges may be used for the evaluation of the low- and the high-energy 
parts. This gauge "arbitrariness" (rather than "invariance") holds only if linear terms in e are 
neglected. This remarkable fact is a cornerstone in the construction of nonrelativistic QED La- 
grangians (NRQED), where most high-energy effective operators are taken from calculations 
carried out in Feynman gauge, but the ultrasoft-scale effects (where the virtual photon energy is 
of the order of the atomic binding energy) are calculated in Coulomb gauge. 



The general structure of the one-loop energy shift E, for non-S states, is 

E = E L + E H = ^ Za)A ^ c2 F{Za). 
7r n 6 



(7) 



The dimensionless function F[Za) depends on the bound-state quantum numbers and has the 
expansion (again, for non-S states) 



F(Za) = A m + (Za) 2 {A 6l In [(Za)~ 2 } + A m ) 



(8) 
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Here, we ignore higher-order terms irrelevant for the current a (Za) 6 -calculation. The indices 
of the constant coefficients A X y correspond to the power X of Za and the power Y of the 
logarithm ln[(Za)~ 2 ]. The function F(Za) is obtained as the sum 



F(Za) = F H (Za, e) + F L (Za, e) , 



(9) 



where F H corresponds to E H and F L corresponds to E L . The sum F(Za) is independent of the 
cutoff parameter e. 

Let us consider the concrete case of the 8D 3 / 2 state. The high-energy part reads 



Fn{8D 3/2 ) 



— + (Za) 2 
20 V ; 



20893 



31 m 



31 



2419200 2520 e 2520 



In I — 

m 



(10) 



The calculation entails techniques familiar from high-energy physics and is described in detail 
in Ref. Il22ll . The term —1/20 is a consequence of the anomalous magnetic moment of the 
electron as described in detail in Chap. 7 of Ref. [23]. The low-energy part is obtained largely 
by atomic-physics calculational techniques and reads 



F L (8D 3/2 ,Za,e) 



^hi£;o(8D) 



31 m 31 

0.024886 + + In 

2520 e 2520 



+ (Za) 2 

where In k is the familiar nonrelativistic Bethe logarithm. The sum is 

F(8D 3/2 , Za, e) =F H (8D 3/2 , Za) + F L (8D 3/2 , Za) 

31 



(Za) 2 m 



(11) 



— --\nk (8D) + (Za) 2 
20 3 V ' K ' 



2520 



In [(Za)' 2 ] + 0.007723 



and thus free of e, as it should be. A comparison of Eqs. © and (fT2l reveals that 

^6o(8D 3/2 ) = +0.007723. 



(12) 



(13) 



The generalization of the term —1/20 — (4/3) lnfc (8D) to an arbitrary hydrogenic state is 
discussed in|XJ Model examples illustrating the cancellation of the cutoff parameter e can be 
found in Appendix A of Ref. Q and in Section 2 of Ref. Q. 



3 Low-Energy Part 
3.1 Orientation 

In this brief review, we will focus on the low-energy part ("ultrasoft" scale, cu ~ (Za) 2 m). 
According to Eq. (3.109) of Ref. [22], this term may be written as 



E T , = -- 



dk k 



k l k? 



a e 



H B -E B + k 



-a e 



-ifc". 



if, 



(14) 
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where the fully relativistic Dirac wave functions ip is used, and the a 1 = 7°7* matrices are 
a noncovariant form of the Dirac 7 matrices. The quantity k = \k\ is virtual photon energy. 
Note that one should actually supplement a principal value prescription in this formula, because 
the lower-lying states generate poles in the integration region k 6 (0, e) with e 2> (Za) 2 m. 
The ensuing problem of the accurate definition of a resonance eigenvalue has been discussed in 
Ref. [26], in the context of two-loop corrections. 

It is well known that the Dirac Hamiltonian H D in Eq. (fT4l) can be transformed to the 
Schrodinger Hamiltonian plus various relativistic corrections, using the Foldy-Wouthuysen 
transformation U, which leads to a transformed Hamiltonian UHdU~ 1 . It is this transforma- 
tion which is the key to the successful identification of all terms which contribute to the low- 
energy contribution in a given order of the Za expansion. Specifically, one may apply that same 



transformation U to the current operator 1 13], 



The result is 



where 



j 3 = a 3 exp(i k ■ r) . (15) 



j 

Uflf- 1 = — + 5f, (16) 
m 



5j 3 ~ £ (1+i (k.r)-l[k ■ r) ) (17) 



m \ V / 2 

1 1 Za j 

— tvP - - — ? —5- {r x ay 



2 m 3 2 m? r 3 

H (k ■ f] ( k x a) — (k x a] . 

2m V A / 2m V J 

Here, we have dropped terms which couple lower and upper components of the Foldy- 
Wouthuysen transformation and lead to vanishing energy shifts within the context of the e- 
expansion. The term proportional to k x a entails a spin flip, but the leading-order current p 3 /m 
necessitates a change in the angular momentum of the electron. Because of angular momentum 
selection rules, we may thus neglect this term in the calculation. 

Alternatively, one might have obtained the current (fTTt by considering a Foldy-Wouthuysen 
transformation of a Dirac Hamiltonian with a vector potential A included, i.e. a transformation 
of a ■ (p — eA) + (3 m + V instead of a ■ p + (3 m + V. The effective current operator then is 
the term that multiplies A in the transformed Hamiltonian. That latter approach is used, e.g., in 
Ref. [27]. The term 

^ f 1+i (k -f) - ~ (k ■ f) 2 ] (18) 



m \ V / 2 



leads to the so-called quadrupole correction. The other current corrections in Eq. (fTTt are con- 
sidered separately and identifies as "relativistic corrections to the current" (see Refs. ifTiiEsll '). 
Of course, these corrections have to be supplemented by relativistic corrections to the Hamilto- 
nian and to the bound-state energy, as well as corrections to the wave function. These corrections 
are well defined and identify all terms relevant in the order a (Za) 6 m. 
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3.2 Dimensionless energy parameter 



One of the advantages of the Foldy-Wouthuysen transformation is that we may carry out all 
further calculations using the nonrelativistic form of the bound-state prop agator. There exists a 



closed-form Sturmian decomposition for all angular components [29, 30, 31], 



ite propag 



\ l / \ l oo T 2l+l I 2r i \ T 2l+l ( 2r 2 



a u \a u J \a u J (k + l) 2 j+i (I + 1 + k - v 

Here, a = l/(Zam) is the Bohr radius, and (k) c is the Pochhammer symbol. The symbols 
L^ +1 denote the associated Laguerre polynomials. The radial integrals can usually be evaluated 
using standard techniques, which leaves the sum over A; as a final problem in the calculation 
(note that k in this context is to be differentiated from the virtual photon energy which is also 
denoted by k in many calculations in this field of research). In|Bj we discuss some important 
properties of typical hypergeometric function encountered which result from an evaluation the 
sums. Convergence acceleration techniques are useful for accurate calculations of relativistic 
effects, because the A;-sum is typically slowly convergent for high virtual photon energies. 

The nonrelativistic Schrodinger-Coulomb Green function for the hydrogen atom [30] reads 



n. 



l 



H 



r 2 )=J29i(ri,r 2 ;v)Y lm (0i^i)Y l * m (6 2 ,<p 2 ). (20) 
i in 



Here, z = E s — k is an energy parameter that involves the Schrodinger energy E s of the 
reference state. When the photon energy k assumes values in the range (0, oo), the argument z 
decreases monotonically from E s < to — oo. The basic idea is now to set t = y^E s /z, so that 
t runs from a value t = 1 for k = to a value t = for k = oo. As a function of Za and the 
principal quantum number n, the quantity t parameterizes the energy argument z as 

. . (Za) 2 m 



, , Za. \ m 
t = t(z) = — J-—. (21b) 
n \ 2z 

Equation (3.126) of Ref. [22] specializes this transformation to the case n = 2. In Refs. (2; 
we encounter the notation 



Eqm 
Q, 



771 

u = nt=(Za) x l-—, (22) 



which enters into Eqs. (fT9l and d20l) . An expansion for large k then corresponds to an expansion 
for small t (see|Bl). 



3.3 Extraction of a finite part 

In many cases, it is necessary, within bound-state calculations, to extract a nonlogarithmic, 
constant term from an integral that contains a variety of divergent contributions in the ultraviolet 
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domain. Here, we discuss possible algorithms which can be used for such calculations, and 
which may eventually be useful for an independent verification of the results reported here. For 
example, the finite parts which lead to the relativistic Bethe logarithms might be alternatively 
extracted using a discretized Schrodinger-Coulomb propagator on a lattice [32]. In this case, 
it is necessary to have a means for extracting nonlogarithmic terms numerically rather than 
analytically. We believe that an independent verification of the nonlogarithmic terms, using 
purely numerical methods, could be a very worthful cross-check of the analytic approach (up to 
the last stage of the calculation where the integrals are evaluated numerically) which was used 
in previous evaluations. 

In order to approach this problem, we first consider a model calculation inspired by the structure 
of the nonrelativistic integrand that leads to the Bethe logarithm in Eq. dTTt below. It consists in 
the evaluation of the integral 



k + h(Za) 2 

= (ZafA - {Zaf h In (^^) (23) 

and in the extraction of the logarithmic and nonlogarithmic terms (e prescription), or in the cal- 
culation of just the nonlogarithmic term (numerical evaluation of generalized Bethe logarithms). 

The first prescription would consist in expanding the integral for large A, and dropping all linear, 
quadratic etc. terms for large A. In this way, we obtain the result 

Let us now consider a variation of the first prescription. We go to atomic units, i.e. scale all 
powers of Za out of the integrand via the transformation k — > (Za) 2 k. We then define A = 
A/(Za) 2 and obtain 

1 ={Za) A [ dkk ' 



o 



--{Zafh 



A - In 



k + h 
h + X 



- {Zafh 



In 



(25) 



where the last form is obtained after dropping the linear and (possibly further) quadratic terms; 
the final result is in full agreement with Eq. d24ll . 

The second prescription is based on the identification A — > e. We first expand the result in 
Eq. (l23l in powers of Za, which entails the replacement 

e_±(Zafh e__ 

(Za) 2 h ~^ (Za) 2 h { ' 
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in the argument of the logarithm in Eq. (1231) . We then expand in e for small e, dropping the linear 
terms. We thus obtain the result 



- (Za > h ln I WfK ) • (27) 



which is equivalent to the result in Eq. (1241) upon the identification A — ► e. This illustrates that 
the expansion in small e after the expansion in Za is actually an expansion for e ^> (Za) 2 m, 
after dropping the linear terms. In particular, all of the above prescriptions lead to the result 

/ ~ (Zafh \n((Za) 2 h) , (28) 

for the nonlogarithmic term. 

Let us now suppose that we can evaluate a general function f(x) only numerically, but that we 
know its expansion for large x, 

f( x )= a n x n + 0(x~ 3/2 ), x-^oo, (29) 

n=— 1,-1/2,0,1/2, ...,rn 

and we wish to evaluate the nonlogarithmic term iV that is generated by the integral 

/ dxf(x)~N (30) 
Jo 

for large A. In the sense of Eq. (l23t . we would have f(x) = x(Za) 2 /[x + h(Za) 2 ] and and 
= (Za) 4 h ln ((Za) 2 h). For the more general case [see Eq. (l2"9tl. we define 

d( x )= a ^ n ( 3l ) 

n=— 1,-1/2, 0, 1 /2,.--: m 



(J/77 Ju 



and 

D( x ) = _! Hx) + ^pr- (32) 

n=—y2,0,y2,...,m 

For arbitrary M, the nonlogarithmic term A^ may then be extracted according to 

N=X l +l 2 +X^ (33) 

with 

-M 



pivi 

lx= dxf(x) , 
Jo 

Z 2 = / da; [f(x) - d(x)} , 

J M 



X 3 = -D(M). (34) 

The sign of the T 3 -term is determined by the necessity to subtract the integral D(x) of the 
subtraction term d(x) at the lower limit of integration x = M. This consideration effectively 
results in three minus signs. Analogous considerations have recently been used in Refs. 13 3l 13411 . 
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4 Numerical Results 



For the states under investigation, the A m coefficients are listed in Tables 1 and 2. For the 
5F 5 /2 state, the result had previously been recorded as 0.002 403 158, and for 5F 7 / 2 , a value 
of 0.008 087 020 had been indicated (see Ref. J20J]). The correction of this result, in the last 
decimal, is beyond current and projected levels of experimental accuracy. For the current brief 
review, we re-evaluate many of the integrals leading to the A 60 coefficients with an enhanced 
number of integration nodes. The two entries in question change by more than the previously 
indicated numerical accuracy. Results for 6F, 7F, 6G, and 7G as reported in Tables 1 and 2 are 
obtained here. In Ref. Ill 811 . we already corrected a computational error for 3Pi/2 as previously 
reported in Eq. (96) of Ref. jH, where a value of -1.14768(1) had been given. As in previous 



calculations (see Refs. I113L 12811 ). certain remaining one-dimensional integrals involving (partial 
derivatives of) hypergeometric functions could only be evaluated numerically. For n = 8, we 
recall the results (H A 60 (8D 3/2 ) = 0.007 723 850 and A 60 (8D 5/2 ) = 0.034 607 492. 

Table 1 : A table of A 60 coefficients for higher excited atomic states with positive 
Dirac angular quantum number k (i.e., j — I — 1/2). All decimal figures shown are 
significant. 

n A m (nP 1/2 ) A eo (nD 3/2 ) A 60 (nF 5/2 ) A 60 (nG 7/2 ) 

2 -0.998 904 402 - 

3 -1.148 189 956 0.005 551 573 

4 -1.195 688 142 0.005 585 985 0.002 326 988 

5 -1.216 224 512 0.006 152 175 0.002 403 151 0.000 814 415 

6 -1.226 702 391 0.006 749 745 0.002 531 636 0.000 827 468 

7 -1.232 715 957 0.007 277 403 0.002 661 311 0.000 857 346 



Table 2: Analog of Table 1 for states with negative k. 

n A 60 (nP 3/2 ) A 60 (nD 5/2 ) A 60 (nF 7/2 ) A 60 (nG 9/2 ) 

2 -0.503 373 465 - 

3 -0.597 569 388 0.027 609 989 

4 -0.630 945 795 0.031 411 862 0.007 074 961 

5 -0.647 013 508 0.033 077 570 0.008 087 015 0.002 412 929 

6 -0.656 154 893 0.033 908 493 0.008 610 109 0.002 748 250 

7 -0.662 027 568 0.034 355 926 0.008 906 989 0.002 941 334 



5 Conclusions 

The challenges of bound-state quantum electrodynamic calculations are associated to the ac- 
curacy of the experimental verifications, to the significance of the theory for the determination 
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of the fundamental constants, and to the conceptual complexity of the calculations which is 
derived from the apparent simplicity of the physical systems under study. The latter aspect is 
developed to full extent only if an accurate understanding is required in higher orders of pertur- 
bation theory. 

Let us briefly discuss possible extensions of this work in addition to the independent verifi- 
cation using a more numerically inspired approach, as outlined in Sec. 13.31 In Refs. 124 13511 . 
calculations of the fine- structure splitting for P states are described which rely on a form-factor 
approach; this would be equivalent to using an effective operator for the high-energy part. Usu- 
ally, NRQED-inspired calculations involve a high-energy part, which takes care of the con- 
tribution of the high-energy virtual photons, and which is given by effective operators, and a 
low-energy part, which is given by photons whose energy is of the order of the electron binding 
energy. The latter integration region is often referred to as the "ultrasoft scale" in the literature 
(see, e.g., Refs. l2lla.l2ll2a.l22D)- The two scales (sometimes three, if one makes an additional 
distinction with regard to electron momenta) require a completely separate treatment and can- 
not be calculated on the same footing. Thus, the introduction of scale- separation parameters is 
required. These cancel at the end of the calculation. For the high-energy effective operators, 
this scale-separation parameter takes the role of an infrared cutoff, whereas for the low-energy 
contributions, the scale- separation parameter gives a natural scale for the failure of the nonrela- 
tivistic ("ultrasoft") approximation to the virtual photons, i.e. it acts as an ultraviolet cutoff. This 
property is characteristic of QED bound-state calculations and is a feature that adds a certain 
twist to the analysis which is not present in usual quantum-field theoretic loop calculations. 



A Mass renormalization and the leading-order result 

In the first articles on the self-energy, the concepts of mass renormalization and covariant inte- 
gration were developed and applied for the first time to the calculation of observable physical 
effects. The leading-order energy shift, of order a(Za) 4 mc 2 , is the sum of a Dirac F\ form- 
factor contribution, an anomalous magnetic moment term (F 2 electron form factor), a vacuum- 
polarization term, and an effect mediated by low-energy virtual photons. The latter contribution 
can be interpreted as an "average virtual excitation energy", as a "logarithmic sum" or, as it 
is most widely called, a "Bethe logarithm" In k . One may therefore point out that the leading 
order-a(Za) 4 -effect summarizes already the core of most properties of the electron mediated by 
the virtual interactions with the quantum fields. The self-energy of a bound electron is actually 
the difference of the self energies of a bound and a free electron, the latter being attributable 
to its mass, wherefore it can be reabsorbed into a redefined physical parameter entering the 
Lagrangian of QED. 

Yet at the same time, it is important to remember that different energy scales enter the prob- 
lem: (a) the atomic binding energy and (b) the relativistic electron rest mass energy scale. The 
separation is already explained in §123 of Ref. [5], and its extension to higher orders in the Za 
expansion is discussed in this brief review. 

The generalization of the a(Za) 4 m-part of the result in Eq. (fT2l to an arbitrary non-S state, 
reads 

„ a (Zafm ( 1 4 , , , ,A 

E ~ - 1 ; - --Infco V . (35) 

7r n 6 \ 2 k (21 + 1) 3 / 



11 



The A i0 -coefficient for a non-S state thus reads (see, e.g., Refs. lEill and lEijl ) 



A. 



■10 



2k (21 + 1) 



- lnA;o(nO, 



(36) 



where k — 2(1 — j)(j + 1/2), with the usual meaning for the bound-state angular momen- 
tum quantum numbers. The Bethe logarithm In ko(nl) is an inherently nonrelativistic quantity, 
whose expression in natural units reads 



In ko(nl) 



n 



2(Za) A m 



P 
m 



(H s - E n ) In 



2 \H S — E n \ 


pi 


(Za) 2 m 


m 



(37) 



Here, <p is the nonrelativistic (Schrodinger) form of the bound-state wave function. 

It took 50 years |EI 0. & HI 01 to advance our understanding from the order a(Za) 4 m e c 2 to 
the order a(Za) 6 m e c 2 . For a successful calculation of the higher-order relativistic corrections, 
it has been essential to master the much more involved analytic, covariant integrations in the 
high-energy part, and to advance our understanding of relativistic and retardation corrections to 
the current. Also, the possibility of handling very involved analytic intermediate expressions by 
computer algebra has become an essential ingredient of the modern calculations. 

Finally, let us remark that for low -lying states, it is possible to evaluate the one-loop effect 
numerically to very high accuracy 114 lLl42l 14311 . even for low nuclear charge numbers. For very 
highly excited states, however, first exploratory work has revealed that the numerical difficulties 
associated with the renormalization are still rather prohibitive in the context of an accurate 
numerical evaluation, because it entails a loss of more than twelve decimals for higher excited 
states. In addition, the wave functions become much more complex, and this inhibits a fast 
convergence of the numerical integrations, leading to the temporary conclusion that analytic 
calculations are still preferable for higher excited states, in the domain of low nuclear charge 
numbers. 



B Asymptotics 

We here discuss the asymptotics of two hyperge ometric functions that are often encountered in 
bound-state calculations of the kind EHElSj considered in the current brief review. The 
first of these functions is 



oo 



$x(n,t)= 2 F 1 (l,-nt,l-nt,£) = -nt^j-^, (38) 

which develops singularities at t = k/n. These correspond to lower-lying states and their 
equidistant spacing is a consequence of the energy parameterization discussed in Sec. 13.21 For 
given n, there are typically n — 1 lower-lying states accessible by dipole decay, and thus n — 1 
singularities at t = m/n, for m — 1, . . .n — 1, as t is in the range (0, 1). The argument £ of the 
hypergeometric typically reads 
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The small-t asymptotics read 

$i(n, t) =1 + nt \n(4t) - [2n + n 2 ((2)] t 2 

+ {n + 4n 2 [1 - ln(4t)] - n 3 C(3) } t 3 



+ 



— n + 4n 2 + 4n 3 C(2) - n 4 ((4) 
3 



t 4 + O [t 5 ln(t)] . 



The other function which is often encountered reads 

$ 2 (n, *, = 2^1(1, -nt, C) 

where 



fc=0 



k — nt 



c 



l-t 
T+t 



Its asymptotics are given by 

$ 2 (n,t) =1 + ntln(2) + 



n 2 C(2)-n 



+ 



+ 



^-2n 2 ln(2) + ^ 3 C(3) 



o o 



t 4 + O [t 5 \n(t)] 



(40) 



(41) 



(42) 



(43) 



Numerical algorithms useful for different t-ranges are discussed in Table 3. For a descrip- 
tion of the combined nonlinear-condensation transformation (CNCT), the reader is referred to 
Refs. flB. 

Table 3: Numerical algorithms used in the calculation of the 
functions $i and $ 2 . 







< t < 0.05 


0.05 < t < 1 


$i(n, 
$2(n, 


*, 
f, 


CNCT 
5 transformation 


power series +recursion 
power series +recursion 
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